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Abstract: In this paper we present a general mathematical construction that allows us to define a para- 



o 

] metric class of i/-sssi stochastic processes (self-similar with stationary increments), which have marginal 
probability density function that evolves in time according to a partial integro-differential equation of 
fractional type. This construction is based on the theory of finite measures on functional spaces. Since 
• the variance evolves in time as a power function, these i?-sssi processes naturally provide models for slow 
^ I and fast anomalous diffusion. Such a class includes, as particular cases, fractional Brownian motion, grey 
' Brownian motion and Brownian motion. 



1 Introduction 



d , The grey noise theory introduced by Schneider (see [14, 15]) leads naturally to a class of self-similar 

^ ' stochastic processes {B/^it), < (3 < 1}. These processes, called grey Brownian motion, provide stochas- 

I— 1| tic models for the slow-anomalous diffusior0 described by the time fractional diffusion equation; i.e. the 
marginal density function of the grey Brownian motion is the fundamental solution of the time frac- 

[ tional diffusion equation (see [16] and [6,7]). This class will be extended to a class {Ba.pit)}, with 

(y-j ' 0<Q!<2, 0</3<l, called "generalized" grey Brownian motion, which includes stochastic models 

\^ ', either for slow and fast-anomalous diffusion. First, we present and motivate the mathematical construc- 

' tion. Then, we show that this class is made up of 7J-sssi processes and contain either Gaussian and 

', non-Gaussian processes (like fractional Brownian motion and grey Brownian motion). Finally, we show 

y—i • how the time evolution of the marginal density function is described by partial integro-differential equa- 

' tions of fractional type. 

i> : 

^P. I We begin introducing some basic concepts and facts. Let X be a vector space over a K-field and 

K*" ■ let {II • lip, p e /} be a countable family of Hilbert-norms defined on it. The space X along with the 
k>( ] Hilbert-norms {| | • | |p, p € /} is said a topological vector space if it carries as natural topology the initial 
■ topologjH of the norms and the vector space operations. We indicate with Xp the completion of X with 
pi ] respect to the norm || • ||p. Let (•, •) denote the natural bilinear pairing between X and its dual space X' . 
We equip X' with the so called weak topology, which is the coarsest topology such that the functional 
{■,x) is continuous for any x £ X. 



Definition 1.1 (Nuclear space). A topological vector space X , with the topology defined by a family of 
Hilbert-norms, is said a nuclear space if for any Hilbert-norm \ \ - Jlp there exists a larger norm \ \-\\q such 
that the inclusion map X„ ^ Xp is an Hilbert- Schmidt operatoi^. 



1 Anomalous diffusion is characterized by the (asymptotic) time power-law behavior of the variance: cr'^(t) ~ f. Namely, 
the diffusion is slow if the exponent 7 is lesser than one, normal if it is equal to one and fast if it is greater than one. 
^The coarsest topology defined on X which makes these functions continuous. 

^An Hilbert-Schmidt operator is a bounded operator A, defined on an Hilbert space H, such that there exists an 
orthonormal basis {eijig/ of H with the property ||j4ei||^ < 00. 
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Nuclear spaces have many of the good properties of the finite dimensional Euclidean spaces M.'^. For 
example, a subset of a nuclear space is compact if and only if is bounded and closed. Moreover, spaces 
whose elements are "smooth" is some sense tend to be nuclear spaces. In the following example we see 
how nuclear spaces could be constructed naturally starting from an Hilbert space and an operator (see 
Kuo [4]). 

Example 1.1. Let H be an Hilbert space and A an operator defined on it. Suppose that there exists an 
orthonormal bases n = 1, 2, . . .} satisfying the following properties: 

1. They are eigenvectors of A; i.e. for any n > 0: Ahn = A„/i„, A„ S R. 

2. {A„}„>o is a non-decreasing sequence such that: 1 < Ai < A2 < • • ■ < A„ 

00 

3. There exists a positive integer a such that: A^" < 00. 

n=l 

For any non-negative rational number p G Q+, we define a sequence of norms {|| • p G Q+} such that: 
\\^\\p = \\APa ^(^H. That is: 

where (•, •) indicates the H inner product. 

Remark 1.1. For any p E Q+, the norm || • ||p is an Hilbert-norm. Indeed, it comes from the scalar 
product: 

00 

{tv)p^Y.^ni^,h„Mhr.). (2) 

71=1 

For any p e Q+ we define: Xp = {£, E H; \\^\\p < 00}. In view of the above remark, Xp is an Hilbert 
space. Moreover, it is easy to see that for any p > q > 0: 

XpCX,. (3) 

We have the following proposition: 

Proposition 1.1. For any p £ Q+, the inclusion map Xp_f_a/2 ^ Xp is an Hilbert- Schmidt operator. 
Proof: we set hf^ = —phn. The collection {/i^, rt = 1, 2, . . .} is an orthonormal bases of Xp. In fact, for 

An 

any positive integers n and to: 

00 00 
— \^\PfhP h,\(hP h.\ 

An An 



00 00 y2p 

(K, Ki)p = X! KiK^hk){hP^, hk) = -Tjrhp-^nAnk = 5nm- 



k=l k=l 



For each ^ G -^^^+0/2; we indicate with z(^) = C £ the inclusion map. Therefore, for any n > 0: 

An 

and thus by hypothesis 



n— 1 n— 1 

Consider the vector space X = Xp. In view of the above proposition X along with the family of 

p6Q+ 

Hilbert-norms {|| ■ ||p , P G Q+} is a nuclear space. 



2 



Let X be a vector space. A continuous map $ : X ^ C is called a characteristic functional on X if 
it's normalized: 

$(0) = 1, 

and positive defined: 

m 

X! Ci^te -?j)Cj > 0, meZ, {c,},= i,...,™eC, {£.i}^=l,...,r,^& X. 

Let X = M". The Bochner theorem [13] states that for any characteristic functional $ defined on R", 
there exists a unique probability measure fi defined on R", such that 

e^(^'«)d//(a;) = $(0, S. G R". 

Let now X be a topological vector space. In the characterization of typical configurations of measures 
on infinite dimensional spaces the so called Minlos theorem plays a very important role. This theorem is 
an infinite dimensional generalization of the Bochner theorem: 

Theorem 1.1 (Minlos theorem). Let X be a nuclear space. For any characteristic functional $ defined 
on X there exists a unique probability measure fi defined on the measurable space {X',B), where B is 
regarded as the Borel a-algebra generated by the weak topology on X' , such that: 



X' 



= eex. (4) 



Characteristic functional on Hilbert spaces can be defined starting from completely monotonic func- 
tionfl In fact we have the following proposition: 

Proposition 1.2. Let F be a completely monotonic function defined on the positive real line. Therefore, 
there exists a unique characteristic functional $, defined on a real separable Hilbert space H , such that: 

This is obvious because completely monotonic functions are associated to non-negative measure defined 
on the positive real line (see Feller [1]). The converse is also true (see Schneider [14,15]). 

2 White noise 

Consider the Schwartz space 5(R). Equip 5(R) with the usual scalar product: 

itv) = / dtatHt), ^,V^S{R). (5) 



We indicate the completion of 5(R) with respect to eq. ^ with 5o(R) = iL^(R). We consider the 
orthonormal system {hn}n>o of the Hermite functions: 

h^{x) = \ H,,{x)e-^''\ (6) 



■*A function F{t) is completely monotone if it is non-negative and possesses derivatives of any order such that: 

dt'= 



(-i)'^^W>o, t>o, fcez+ = {o,i,2,...} 
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where Hn{x) = (— l)"e^ {d/dx)"e ^ are the Hermite polynomials of degree n. Let A be the "harmonic 
oscillator" operator: 

A^-^+x^ + l; (7) 

A is densely defined on 5o(K) and the Hermite functions are eigenfunctions of A: 

Ahn = A„/i„ = (2n + 2)hn, n = 0, 1, . . . . 

We observe that 1 < < Ai < . . . < A„ and A^^ < oo. We are in the condition of Example 11.11 
Therefore, for any non-negative integer p, we can define: 

/ oo \ 1/2 

\n=0 / 

where 1 1 • 1 1 indicates the L"^ norm. The Schwartz space 5(]R) could be then "reconstructed" as the projective 
limit of the Hilbert spaces 5p(M) = G M\p < oo}. That is: 

S{R) = fl Sp{R). (8) 

p>0 

Therefore, the topological Schwartz space, with the topology defined by the || • ||p norms, is a nuclear 
space. Since 5(R) is a nuclear space, we can apply the Minlos theorem in order to define probability 
measures on its dual space iS'(M). Consider the positive function F{t) ~ e~*, t > 0. It is obvious that 
F is a completely monotone function. Therefore, the functional $(^) — -^dl^lP), C G •^^(IR); defines 
a characteristic functional on iS(M). By Minlos theorem, there exists a unique probability measure /i, 
defined on (5'(R),S), such that: 

e''<"'«>d/i(tj) = e-ll«ll', ee5(K). (9) 

5'(R) 

The probability space (iS'(R), S, /i) is called white noise space and the measure ^ is called white noise 
measure, or standard Gaussian measure, on 5'(R). 

Consider the generalized stochastic process X, defined on the white noise space, such that for each 
test function (p G 5(M): 

X{ip){-)^{;^). (10) 

Clearly, for any ip e S(M.}, X{(p) is a Gaussian random variable with zero mean and variance E{X{ip)'^) — 
2||v3|p. Moreover, for any (/?, G S{R): 

E{X{ip)Xi^)) = 2{ip,(P), (11) 

where E{w) indicates the expectation value of the random variable w. We refer to the generalized process 
X as the canonical noise of (5'(M), S, /i). 

Remark 2.1. In view of the above properties the process X is a white noise [4], and this also motivate 
the name "white noise space" for the probability space {S'{R),B,^). 

We have the following: 

Proposition 2.1. For any h e iL^(]R), X{h) is defined almost everywhere on iS'(M). Moreover, it is 
Gaussian with zero mean and variance 21 
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Proof: we indicate with (iL^) = £2(5'(R),^). Clearly, for any ^ € we have that X(^) e (il^) and: 

um\!^^) = Eixio') = m\\h- (12) 

For each /i € 'f^^(R), there exists a sequence {^n}nGN of 5(R)-elements which converges to h in the iL^(R)- 
norm. Therefore, from eq. (fT2l) . the sequence {-'f (CTi)}riGN is Cauchy in (iL^) and converges to a limit 
function X{h), defined on S'{R). □ 

The latter proposition states that for every sequence {/t}tgR of £-^(M)-functions, depending continu- 
ously on a real parameter i e K, there exists a Gaussian stochastic process 

{Y{t)}teM = {X{ft)}tm, (13) 

defined on the probability space (5'(M), S, /i), which has zero mean, variance E(Yt)'^ = 2||/t|p and 
covariance E{Yih)Y{t2)) - 2{ft,Jt,). 

Remark 2.2. Observe that if W{x), a; G M, is a Wiener process defined on the probability space {fl, J-, P), 
then the functional 

X{ip)= [ ip{x)dW{x), (peii^iR), (14) 



is a white noise on the space {n,J^,P). Therefore, if we indicate with l[o,t)(a;), t > 0, the indicator 
function of the interval [0,t), the process 



is a one-sided Brownian motion. 



X(l[o,t))= / dWix)=Wit), t>0, (15) 



Example 2.1 (Brownian motion). Let X be a white noise defined canonically on the white noise space 
(5'(R), S, /i). Looking at eq. ifTS]) . its natural to think that the stochastic process 

{B(i)}t>o-{X(l[o,t))}t>o, (16) 

defines a "standard" Brownian motior@. Indeed, the process {X(l[o, t))}t>05 is Gaussian with covariance: 

E [X(l[o,t))X(l[o,.))] = 2(l[o,t), l[o,s)) = 2min(t,s), t,s>0. 

Example 2.2 (Fractional Brownian motion). The stochastic process: 

{B^/2{t)}t>o = {XiUt)}t>o, 0<a<2, (17) 

where 

/o,t(x) = ^^((i-x)7^-(-x)7^), x+=max(x,0), (18) 

and 

(r(a + l)sin^)^ 

is a "standard" fractional Brownian motion of order H ~ a/2 (see Taqqu [18]). 



^With the word "standard" Brownian motion we mean that E{B{1) ) = 2. 
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3 Grey noises 



We have seen that white noise is a generalized stochastic process X defined canonically on the white noise 
space (5'(R),S, /i), with space of test functions We have remarked that the white noise could also 

be defined starting from stochastic integrals with respect to the Brownian motion. In this case the space 
of test function turns out to be the space of integrands of the stochastic integral. Then, the Brownian 
motion B{t) could be obtained from the white noise by setting B{t) ~ X(l[o.t)). We want to generaHze 
the previous construction in order to define a general class of 77-sssi processes which includes, Brownian 
motion, fractional Brownian motion and more general processes. 

Consider a one-sided fractional Brownian motion {Ba/2{t)}t>o with self-similarity parameter H = a/2 
and < a < 2, defined on a certain probability space (fJ,JF, P). The fractional Brownian motion has a 
spectral representation [18]: 

Bo^Mt) = VC{a) f -^^——^\x\^dB{x), t>0, (20) 



wheTedB{x) is a complex Gaussian measure such that c?i3(a;) = dBi{x)+idB2{x) with dBi[x) — dBi{—x), 
dB2{x) = —dB2{—x) and where Bi and B2 are independent Brownian motion. Moreover, 

C(a) =r(a-Hl)sin^. (21) 

We observe that 

1 e''^ — 1 ~ 

where we have indicated with f{x) the Fourier transform of the function / evaluated on x G M: 

J{x) = J(/) (x)^^ f e^^yf{y)dy. (23) 



In view of eq. l(22|) we have: 



Ba./2{t) = VC{a) / l[o,t)(x)|x|^dB(x). (24) 



Therefore, if one defines a generaHzed stochastic process X such that for a suitable choice of a test 
function ip 



X^iip) = y/C{^ / ^{x)\x\^dBix), (25) 

JR 

one can write: 

S„/2(0 =^«(l[o,t)), t>0. (26) 
Remark 3.1. The space of test function can be the space 

A. = {/ e £2(R); ml ^ C{a) f \f{x)\'\x\'-"dx < 00}, (27) 

which coincides with a space of deterministic integrands for fractional Brownian motion (see Pipiras and 
Taqqu [10,11]). 
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Consider now the Schwartz space 5(R) equipped with the scalar product: 



{^,v)a=C{a) U{x)?l{x)\x\^~''dx, ^,77e5(R), 0<a<2, (28) 



where C{a) is given by eq. ((2T|). This scalar product generate the a-norm in eq. l(27|) . We indicate with 
5^"^(M) the completion of 5(R) with respect to eq. ([28l) . 

Remark 3.2. If we set a = 1 in eq. ((28l). we have C(l) = 1 and: 



{^,v)i = / axMx)dx = / ay)v{y)dy, (29) 

SO that, we recover the £^(R)-inner product. Moreover, 5q^''(]R) = Sa{M.) = iL^(R). 

Starting from the Hilbert space (Sq"\m.), \ \ ■ ||„), it is possible to reproduce the construction of Example 
11.11 Then, the space turns out to be a nuclear space with respect to the topology generated by 

the a-norm || • \\a and an operator Here we just say that the main ingredient are the Generalized 

Laguerre polynomials: 

^;i(x) = ^|;(e-^x-^), ,>-l,x>0, (30) 

where n is a non-negative integer. They are orthogonal with respect to the weighting function x"'e~^, 

x^e-^Llix)LZ^{x)dx^ W in (31) 

1 (n + 1) 

and satisfy the Laguerre equation: 

'=d^ + ^^+^-''^d^)^ni^^^-''^'r.ix^- (32) 

Using eq. ((3T1) . it is easy to show that the sequence of functions {/i"}nGZ+ defined by: 



'■2n+l 

is an orthonormal bases of 5g"' (R) with the choice: 



^2n+i(2;) = &n,Qe ^'^/^xLI "^^{x^), neZ+, 



(33) 



(34) 



C{a)r{n + l-a/2)J ' \C{a)T{n + 2 ~ a/2) 

Then, using eq. l[32|) . one can show that the orthonormal bases {h^}nez+ is a set of eigenfunction of an 
operator A^°'\ defined on Sq"\m), with eigenvalues aI"-* = 2n + 2 — a + 1. 

Remark 3.3. We recall the well known relationships between Laguerre and Hermite polynomials: 

H2nix) = (-l)"22»n!L,T'/'(x2) 

(35) 

i?2„+i(.x) = (-l)"22"+in!xLy'(x2). 

In view of the above relations, when a = 1 the orthonormal bases {h"}niEZ+ reduces to the Hermite bases 
of iL^(R) eq. ([6]), which is preserved under Fourier transformation. 
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By Proposition [121 starting from a completely monotonic function F, we can define characteristic func- 
tionals on 5(M) by setting = -P'(IICIIa)- Then, we could use Minlos theorem in order to define 

probability measures on 5'(M). We consider the real valued Mittag-Leffier function of order /? > 0: 

W-Ej.^^, xeM. (36) 

n=0 ' 

It is known that the function Fj^ii) = Ep{—t), t > 0, is a completely monotonic function if < /3 < 1 [8]. 

For example if /? = 1 we recover Fi{t) = e"*. Therefore, the functional $a,/3(?) = Fp{\\^\\l), ^ G 5^"^(M), 
defines a characteristic functional on iS(]R). By Minlos theorem, there exists a unique probability measure 
^J,a,f3, defined on {S'{M.),B), such that: 



e 



^■^<'-'^U^,^A^) ^ F(,m\l), Ce5(]R). (37) 



When a — P and < f3 < 1, the probability space {S' (R), B, fip.p) is called grey noise space and the 
measure ^i3,p is called grey noise measure (see Schneider [14,15]). In this paper, we focus on the more 
general case < a < 2 and we call the space {S'{R),B, fia,p) "generalized" grey noise space and fia,0 
"generalized" grey noise measure. 

Definition 3.1. The generalized stochastic process Xa^p, defined canonically on the "generalized" grey 
noise space {S'{M.),B, ^ia,i3), is called "generalized" grey noise. Therefore, for each test function ip e 5(R).- 

Xo.Av){-)^{;V)- (38) 
Remark 3.4. By the definition of "generalized" grey noise measure eq. ((37)) . for any e 5(R), we have: 

ij(e'^y^...M)^£;^(_y2||^||2)^ yGM. (39) 

Using eq. ([39l) and eq. ([36]) it easy to show that the "generalized" grey noise has moments of any order: 



la ■) 



for any integer n > and ^ € 5(IR.). It is possible to extend the space of test functions to the whole 
5^"^(R). In fact, for any ^ € S{M) we have X^AC) ^ ('^^) = '^'^(^'(M), /i„,/3). Thus, for any h € 5^"^(M), 
the function XaA^) defined as a limit of a sequence Xa^p{S^n)i where {Cn} belong to 5(R). Therefore 
we have the following: 

Proposition 3.1. For any h e 5q"''(M), Xa i3{h) is defined almost everywhere on iS'(M) and belongs to 

Summarizing: the "generalized" grey noise is defined canonically on the grey noise space [S'iR),B, fia,f3) 
with the following properties: 

1. for any h e 5q"''(M), Xa,p{h) is well defined and belong to (il^). 

2. Eie'y^".'^^''')) = Efj{-y^\\h\\l) for any yeR. 

3. E{X^.p{h)) = and EiX^Ahf) = Y\^^)\Ml 
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4. For any h and g which belong to 5q"''(M), one has 



1 



If we put /3 = 1, the measure := is a Gaussian measure and Xa^i := Xa is a Gaussian noise. 
In fact, for any h e 5g"''(R), the random variable Xa{h) is Gaussian with zero mean and variance 
E[Xa{hY) = 2||ft,||^ (see eg. [39|) . When a = 1, Xa reduces to a "standard" white noise (see Remark 13.21 
and Remark [33| . Moreover, for any sequence {/t}teK of 5g"^ (M)-functions, depending continuously on 
a real parameter t g R, the stochastic process Y{t) ~ Xa{ft) is Gaussian with auto-covariance given by 

eq. Igl]) 

E{Y{t)Y{s)) = E{Xa){ft)XMs) = {ftjs)c. + {ftjs)a- (42) 

Example 3.1 (Fractional Brownian motion). For any t > the function l[o,() belongs to 5g"''(R). In 
fact, it is easy to show that ||l[o,t)||Q < oo when < a < 2 and 

l|l[o.)ll^ = ^//^^(l-costx) = t". (43) 

Therefore, we can define the process: 

^a/2W =^a(l[0,t)), t>0. (44) 

The process B^/^ii) is a "standard" fractional Brownian motion with parameter H ~ a/2. Indeed, it is 
Gaussian with variance E{Ba/2{t)'^) = 2||l[o,t)lla = and auto-covariance: 



E{Ba/2it)Ba/2is)) = (1 [0,t) , 1 [0,s) )q + (l[0,t), l[0,s))^ 

f dx- — ; — — (1 — costx + 1 — cos — 1 + cos(t — s)a;) 
= e + s°'-\t-s\" =ja{t,s), t,s>0, 



2-K 



which is the fractional Brownian motion auto-covariance. 

In view of the above example, Xa could be regarded as a fractional Gaussian noise defined on the space 
(5'(R),S,Ma). 

Example 3.2 (Deconvolution of Brownian motion). The stochastic process 

{B{t)}t>a = {Xa{go.,t)}t>o, (45) 
where, for each t > 0, the function ga,t is defined by: 

9o,A^) = —j=^\Q,t){x){ix)^ , (46) 

is a "standard" Brownian motion. Indeed, it is Gaussian, with zero mean, variance 

E{B{tf) = 2 / |x|i-"|T[o,t)(x)na;r-ida; = 2 / |T[o,t) (a;)pdx = 2i, (47) 

and autocovariance: 

E{B{t)B{s)) = / (^){x)\o^,){x) +%~i{x)\o,t){x)^ dx = 2min(t, s). (48) 
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Remark 3.5. The representation of Brownian motion in terms of the fractional Gaussian noise eq. I|45p 
corresponds to a particular case of the so called deconvolution formula, which expresses the Brownian 
motion as a stochastic integral with respect to a fractional Brownian motion of order H = a/2 (see [12]). 
More generally, we can represent a fractional Brownian motion B^/2{t) of order i7 = 7/2,0<7<2in 
terms of a fractional Gaussian noise of order a, which corresponds to a representation of in terms 
of a stochastic integral of a fractional Brownian motion -8^/2 of order H = a/2, < a < 2 (see example 
below). 

Example 3.3 (Deconvolution of fractional Brownian motion). The stochastic process, 

{B^/2it)}t>0 = {^a(5a,7,t)}t>0, (49) 

where: 

lch)~ o--, 
= y ^^l[o,t)(a;)(ix) 2 , 0<7<2, (50) 

is a "standard" fractional Brownian motion of order H = j/2. 

We consider now the general case 0<Q!<2,0</3<1. 
Definition 3.2. The stochastic process 

{Ba,pit)}t>0 = {-'^Q,;3(l[0,t))}t>0, (51) 

is called "generalized" (standard) grey Brownian motion. 

The "generalized" grey Brownian motion Ba^p has the following properties which come directly from the 
grey noise properties and eq. (^3]) : 

1- Ba^piO) ~ almost surely. Moreover, for each t > 0, E{Ba.i3{t)) ~ and 

2. The auto-covariance function is: 

E{B^,p{t)B^,p{s)) = 7a,/3(i, s) = Y\fTT) + - 1^ - ■ (53) 

3. For any t,s >0, the characteristic function of the increments is: 

£;(e*«(^°.f(*)~-«°.'3(-^») =£;/3(-y2|t-s|"), y € R. (54) 

The third property follows from the linearity of the grey noise definition. In fact, suppose < s < i, we 
have y{Ba,0{t) - B^As)) = yXa,p{l[o,t) - l[o,s)) = ^a,/3(yl[s,t)), and ||2/l[s,t)||a = y^{t- s)". All these 
properties are enclosed in the following: 

Proposition 3.2. For any < a < 2 and < /? < 1, the process Ba,/j{t), t > 0, is a self-similar with 
stationary increments process (H-sssi), with H — a/2. 

Proof: This result is actually a consequence of the linearity of the noise definition. Given a sequence of 
real numbers {6'j}j=i,...,n, we have to show that for any < < i2 < ■ • ■ < and a > 0: 

E exp{iJ2^jB^Aat,)) = ^ exp(z ^ 0,a^B„,^(<j)) 
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The linearity of the grey noise definition allows to write the above equality as: 



E 



exp iXQ,/3(^6'jl[o, 



atj)) 



= E 



cxp I iXa,f}{a'i ^6'jl[o,t,)) 



Using eq. I|39p we have 



which, because the complete monotonicity, reduces to 



In view of the definition eq. ([281) and eq. l(22|) . the above equality is checked after a simple change of 
variable in the integration. In the same way we can prove the stationarity of the increments. We have to 
show that for any ft, G M: 



E 



cxp iJ2^ABa.p(tj + h) - BaA^)) 



= E 



cxp i^OjiBaAt]) 
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We use the linearity property to write: 



E 



exp iXa,i3 



'j^lh,tj+h)) 



= E 



exp iXa^pi 



By using the definition and the complete monotonicity, we have: 







jl[/l,tj+h)||Q 



3 



which is true because: 



\h,tj+h){x) 



1 e'- 



27r ix 



□ 



In view of Proposition 13. 2[ {Bajj{t)} forms a class of i?-sssi stochastic processes indexed by two 
parameters < a < 2 and < /? < 1. This class includes fractional Brownian motion (/3 = 1), grey 
Brownian motion {a = (3) and Brownian motion {a = (3 = 1). In Figure [llwe present a diagram which 
allows us to identify the elements of the class. The long-range dependence! domain corresponds to the 
region 1 < a < 2. The horizontal fine represents the processes with purely random increments, that is, 
processes which possess uncorrelated increments. The fractional Brownian motion is identified by the 
vertical line (/3 = 1). The lower diagonal line represents the grey Brownian motion. 



4 Master equation and concluding remarks 

The following proposition characterizes the marginal density function of the process {Ba,fj{t), i > 0}: 

Proposition 4.1. The marginal probability density function fa,(3{x,t) of the process {Ba^pit), t > 0} is 
the fundamental solution of the "stretched" time-fractional diffusion equation: 

u{x,t)=uo{x) + ^ ^s'^/l'-^ (^t9 - '^u{x,s)ds, t>Q. (55) 

Proof: eq. (|54| (with s = 0) states that fa,i3iy,t) = £^/3(— y^t"). Using eq. ((SB)) . we can show that the 
Mittag-Leffier function satisfies 

where we have used the change of variables s' = s"'/^ . Thus, faj{x,t) solves eq. l(55|) with initial condi- 
tion uo{x) = /q^^(x,0) = 5{x). □ 

We refer to eq. ([SS)) as the master equation of the marginal density function of the "generalized" 
grey Brownian motion. Therefore, the diagram in Figure [1] can be also read in terms of partial integro- 
differential equation of fractional type. When a = f3 and < /3 < 1, we recover the time-fractional 
diffusion equation of order P (lower diagonal line). When /3 = 1 and < a < 2, we have the equation 

^An H-sssi process is said to possess long-range dependence if the discrete process of the increments exhibits long-range 
dependence. That is, if the increments autocorrelation function tends to zero like a power function and such that it doesn't 
result integrable [18]. 
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of the fractional Brownian motion marginal density, that is the equation of a stretched Gaussian density 
(vertical line). Finally, when a = /3 = 1 we find the standard diffusion equation. 



From Proposition 14.11 it follows that the parametric class {Ba,i3{t)} provides stochastic models for 
anomalous diffusions described by eq. l(55|) . Looking at eq. l(52|) and eq. l(53|) . which describe the variance 
and the covariance function respectively, it follows that: 

When < a < 1, the diffusion is slow. The increments of the process Ba,f3{t) turn out to be negatively 
correlated. This implies that the trajectories are very "zigzaging" (antipersistent). The increments form 
a stationary process which does not exhibit long-range dependence. 

When a ~ 1, the diffusion is normal. The increments of the process are uncorrelated. The trajectories 
are said to be "chaotic". 

When 1 < a < 2, the diffusion is fast. The increments of the process Ba,p{t) are positively correlated. 
So that, the trajectories are more regular (persistent). In this case the increments exhibits long-range 
dependence [18]. 



The stochastic processes considered so far, governed by the master equation l(55|) . are of course Non- 
Markovian. We observe that non-Markovian equations like eq. (|55|) are often associated to subordinated 
stochastic processes D{t) = B(l{t)), where the parent Markov process B{t) is a "standard" Brownian 
motion and the random time process l{t) is a self-similar of order H = (3 non-negative non-decreasing 
non-Markovian process. For example, in Kolsrud [5] the random time l{t) is taken to be related to the 
local time of a d = 2(1 — /3)-dimesional fractional Bessel process, while in Meerschaert et al. [9] (see 
also Gorenfio et al. [2] and Stanislavsky [17]), in the context of Continuous Time Random Walk, it 
is interpreted as the inverse process of the totally skewed strictly /3-stable process. Heuristically, our 
stochastic process {Ba^pit), i > 0} cannot be a subordinated process (for example if /3 = 1 it reduces to 
a fractional Brownian motion). Therefore, here we provided an example of a class of stochastic models 
associated to time-fractional diffusion equations Hke eq. l(55|) . which are not subordinated processes. 

It is important to remark that, starting from a master equation which describes the dynamic evolution 
of a probability density function /(x,t), it is always possible to define an equivalence class of stochastic 
processes with the same marginal density function f{x,t). All these processes provide suitable stochastic 
models for the starting equation. In this paper we focused on a subclass {i?Q.,/3(i): * ^ 0} associated to 
the non-Markovian equation eq. ((55| . This subclass is made up of processes with stationary increments. 
In this case, the memory effects are enclosed in the typical dependence structure of a fl"-sssi process 
eq. l(53|) : while, for instance in the case of a subordinated process, these are due to the non-Markovian 
property of the random time process. 

It is also interesting to observe that the "generalized" grey Brownian motion turns out to be a direct 
generalization of a Gaussian process. Indeed, it includes the fractional Brownian motion as particular 
case when (5 = 1. Moreover, for any sequence of real numbers {6'i}i=i, ...,„, if one considers the collection 
{Ba,p{ti)^ . ■ . , Ba.pitn)} with < < f2 <•••<*«, it is easy to show that: 

E fcxp(i^0jB„,^(tj)) j =E fexp M^a,/3 ( X] '^jl[o,t,) 



Ep I -r(/3+l)i^0,0,7a,/3(i»,t,) I , (56) 

where ^a,i3 is the the autocovariance matrix eq. ((53|) . It is clear that, fixed /3, the "generalized" grey 
Brownian motion is defined only by its covariance structure. In other words, Ba^p{t) provides an example 
of a stochastic process characterized only by the first and second moments, which is a property of Gaussian 
processes. 
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